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Abstract 

Non-equilibrium phase transitions of a scalar field in an expanding 
spacetime are discussed. These transitions are shown to lead, for 
appropriate potential energy functions, to a biased choice of vacuum 
structure which can be analytically described using percolation theory. 
The initial domain wall networks, which form between different vacua, 
are evolved in time by computer simulation and their behavior in time 
analyzed. It is shown that, unlike systems in thermal equilibrium, 
domain walls formed in biased systems persist for only a short time 
before decaying exponentially away. This result opens the door to a 
complete re-analysis of domain walls in cosmology. 



1 The "Structuron" Field 

We begin by considering a real scalar field, </>, which lives in an expanding 
Friedman- Robinson- Walker (FRW) spacetime. We will assume that this field 
interacts with other fields by gravitation alone, and interacts with itself via 
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a potential energy of the form 



where 



V(</>) = V + 1 



~2 « H , 
V < v 4 



and Hi is the initial Hubble parameter during the epoch of inflation. One 
minimum of this potential occurs at <f> = ttv. Expanding V around this 
minimum gives 

Defining mass m = and coupling A = then the conditions in equa- 
tion (2) imply that 

m<^Hi, (4) 
A < 1. 

It follows that field <f> is not in thermal equilibrium. Many such fields appear 
naturally in various "beyond the standard model" scenarios. Here we will 
simply discuss such a field in the abstract, referring to it by the generic 
name "structuron" since, as we will see, it can play a role in creating large 
scale extra-galactic structure. The behavior of a structuron is determined 
by considering the quantum mechanics of a real scalar quantum field, </>, in 
expanding FRW spacetime. This quantum field can be expanded as 

I = <M + 4> q . (5) 

Field <f> c satisfies the classical equation of motion 

dV 

<f> c + 3i^ c + — = 0. (6) 
o<p c 

During inflation, when H = H{ } and long afterward during the expansion 
epoch, the conditions in equations (4) imply that one can drop the potential 
energy term to leading order. The result is that, to leading order 



V (7) 
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where V is a completely arbitrary constant. To the next order, the classical 
field develops a tiny damped velocity 



which can safely be ignored. The quantum operator, <f) q , and its physical ef- 
fect, is a more complicated issue. It is well known that the quantum fluctua- 
tions associated this field lead, during the inflationary epoch, to the formation 
of a "weakly inhomogeneous, quasi-classical" scalar field, $. This consists, 
just at the end of inflation, of a quantum induced zero mode, plus a part 
which can be Fourier decomposed into wavelengths satisfying the condition 

H- 1 < X < L. (9) 

Here H~ x is the radius of the event horizon at the end of inflation and L 
is the radius of the Universe. During the expansion phase, the radius of 
the event horizon increases with time. At redshift z, this radius is given by 
! c (z) = H^z) -1 . As wavelengths of the quasi-classical field come within the 
growing horizon, they rapidly decay in amplitude and can be ignored. At 
any given redshift z, only those wavelengths satisfying 

l c (z) < X < L (10) 

remain "frozen" and compose the quasi-classical field. Field $, at any redshift 
z, has an interesting spatial decomposition that will be useful in our later 
analysis. To discuss this we first must introduce the two-point correlation 
function, £(/), defined by 

= (o\Mx + i)Ux)\o) (ii) 

where the wavelengths in the quantum operators are to be cut off as in equa- 
tion (10) and / is arbitrary. The translational invariance of the Bunch-Davies 
vacuum |0) implies that £ is independent of coordinate x. Expression (11) 
can be evaluated, and for l c <C I <C L } we find that 
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This result is completely dominated by the long wavelength modes satisfying 
/ < A < L. Secondly, we must introduce a measure of the random fluctuations 
between two points, A(7), which is defined by 

*<'>=IH- (13) 

This function is completely dominated by the short wavelength modes satis- 
fying l c (z) < A < /. Note that 

«') + A(/) = A ln ^ (14) 

which is the mean square value of all the fluctuations in $. We are now in 
a position to decompose field $, as alluded to above. Consider the universe 
at some fixed redshift z. Let us "coarse-grain" the universe into spheres of 
radius /. To first approximation, all of these spheres are correlated by the 
longer wavelength modes, which they share in common. Therefore, they all 
share an approximate zero mode background, 23/ , defined by 



Si = +V*(0- (15) 

However, within each sphere there are fluctuations due to the shorter wave- 
length modes. Since these fluctuations arise in a causally independent way 
from one sphere to the next, each sphere exhibits an independent fluctuation, 
T\, whose value is determined probabilistically from the Gaussian distribu- 
tion function 

P(Ti) = . 1 e "A). (16) 
Note that the root mean square value of this distribution is given by 



rms 



^) = + VA(Z). (17) 



Unlike B\ which is a spatial constant, T\ varies randomly over the coarse- 
grained manifold. It follows that for any length scale /, the weakly inhomo- 
geneous, quasi-classical scalar field can be decomposed as 

$(x) = B, + F,(x). (18) 
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Figure 1: Distribution function, P(^Fi) and potential, V{cf>). 

At any point in space, the relationship between the total shared zero mode 
background V + 23/, the quantum fluctuation T\ and the potential energy 
is shown in Figure f. As long as the potential energy in equation (6) can 
be ignored, no phase transition takes place. However, after a long time, the 
redshift is such that the kinetic energy density becomes comparable, and then 
smaller than, the potential energy. Let us say that they become comparable, 
and the potential energy no longer ignorable, at redshift z t . At this time, at 
each point in space, the held must roll either toward vacuum (+) or toward 
vacuum ( — ). It is clear from Figure I that the the probability that it roles 
to the (+) vacuum is given by 

p= / dTxPiTx). (19) 
Jo 

It is important to note that, since classical background V is arbitrary, p can 
take any value in the range < p < 1. In general, p ^ 1/2. Similarly, 
the probability that the held rolls to the ( — ) vacuum is given by I — p. 
Putting everything together, we see that at the time of the phase transition, 
z t} in any spatial sphere of radius / (that is, at any coarse-grained point) 
the system settles into vacuum (+) with arbitrary probability p } or vacuum 
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( — ) with probability 1 — p } and that the choice of the vacuum in any two 
spheres is statistically independent [1]. Now between any two neighboring 
spheres which have a different vacuum choice, the scalar field <f> will form a 
topological kink; that is, a domain wall. Hence, at z t an initial domain wall 
network will form over all space. Can we determine the spatial structure of 
this initial domain wall network? Yes, using percolation theory. 

2 Percolation Theory 

Percolation theory predicts statistical characteristics and topological prop- 
erties of a "typical" vacuum pattern formed on a lattice. If we partition 
space into a cubic lattice, with each cube having volume / 3 , then percolation 
theory becomes applicable to the prediction of the initial distribution of do- 
main walls. In the following, we will always restrict our discussion to three 
dimensional, cubic lattices. Percolation theory tells us that there exists some 
critical probability, p c . If the population probability for a vacuum exceeds 
this critical value, that vacuum will percolate the lattice; that is, one can 
trace a path from one face of the lattice to another without crossing a do- 
main wall. If the bias probability of a vacuum is less than the critical value, 
however, that vacuum will not percolate, and domain wall "bags", or clus- 
ters, of this vacuum will form. Whether one, both or neither vacua percolate 
depends on the relative values of p, 1 — p and p c . 

It has been shown that p c = 0.31 1 is the critical probability for a cubic 
lattice in three dimensions [2]. Thus, for p < p C} the (+) vacua are in finite 
clusters while the ( — ) vacua sites lie predominantly in a large percolating 
cluster, since necessarily 1 — p > p c . It follows that the associated domain 
walls are relatively small, forming around the compact boundaries of the 
finite clusters. Figure 2 shows this behavior in the case of a small bias 
probability, p = 0.1 < p c . Note that there exist only finite, disconnected 
domain wall bags, the number and size of which can be shown to grow as 
p approaches p c from below. For the case that both p and 1 — p exceed 
p C} both vacua will percolate the lattice, and infinite (that is, lattice sized) 
domain walls separating the vacua will form. A small number of domain wall 
bags, disconnected from the percolating cluster, also form but their number 
decreases as p approaches 0.5. Figure 3 shows this "infinite" domain wall 
structure in the limiting case of p = 0.5. The crucial lesson to be learned 
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Figure 2: Initial distribution of domain walls in three dimensions with bias 
probability, p = 0.1. Here we show a (10A) 3 grid. 




Figure 3: Initial distribution of domain walls in three dimensions with bias 
probability, p = 0.5. For clarity we show here only a (5A) 3 grid. 



7 



from all this is that when both vacua percolate, the topology of the post- 
transition vacuum is one of long, convoluted domain walls stretching across 
all space, whereas, when p < p C} the vacuum is composed of small, compact 
domain wall bags. 

Percolation theory allows one to give a reasonably accurate mathematical 
description of the number of finite s-clusters, their radius and the size of their 
boundary. Here s denotes the number of neighboring lattice sites that are 
occupied by the same vacuum. Let n s be the probability that a given lattice 
site is an element of an s-cluster. This is a fundamental quantity, given 
by the ratio of the total number of s-clusters, N s , over the total number of 
lattice sites, N. An analytical expression for this quantity has been found 
using scaling arguments and Monte Carlo simulations [3]. The result is 



.05013^ e --6299(^K[(^K + 1.6679] 5 (2Q) 



where r = 2.17 and a = .48. Similarly, the average radius of gyration for an 
s-cluster R s (p) } for p < p c and s > s^, is found to be 

R s (p) = .702(p c - P y 322 s- 55 A, (21) 

where 

/ .3ii y- 08 , x 

It can also be shown, for p < p c and s>5, that every s-cluster has a boundary 
composed of 

t s = (±j^J s (23) 

sites. One can easily check using the formula for n s that, on a given lattice, 
the number of s-clusters falls rather quickly with growing s. Hence there 
is an s max such that the total number of s ma:E -clusters is 1. In other words, 
formation of clusters with s much larger than s max is extremely improbable. 
This means that on a given lattice there exists an upper statistical cut-off on 
the size of observable clusters.. 

Using the above formulae, as well as other results found in [3], it is pos- 
sible, for a fixed lattice and a given value of p } to compute the surface area 
for the domain walls between the (+) and the ( — ) vacua. In particular, we 
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have attempted to compute the domain wall surface area for the three di- 
mensional cubic lattice case. However, this calculation is very difficult when 
there are two percolating vacua and hence, we are restricted to a calculation 
for p < p c . This difficulty arises from the fact that the associated large do- 
main walls are highly fractalized and hence their surface area is difficult to 
characterize. The finite size of the lattice further reduces the value of p for 
which we can perform a reasonable calculation to p<0.25. Furthermore, for- 
mula (20) for n s is not very accurate for p < 0.175. This puts a lower bound 
on the calculation. A very reasonable calculation can be performed in the 
range 0.175 < p < 0.25 but, even in this range, we estimate an error of about 
10%. Our results for the domain wall surface plotted, along with the 

10% error bars, as curve (a) in Figure 4. These results can be checked, and 
extended to any values of p } simply by letting the computer evaluate wall 
surface area. The solid line (b) in Figure 4 represents the mean initial surface 
area generated for < p < 0.5. The initial area grows monotonically with 
p, approaching the maximal value of 1.5 X V at p = 0.5. Note that in the 
region where we can compare the percolation prediction with the computer 
experiment, there is good agreement. 

We conclude that, by using percolation theory, we can analytically de- 
termine the structure of the initial network of domain walls. The principal 
result is that for p < p c the domain walls form small, compact bags whereas 
when p > p c the domain walls become very long and stretch across the entire 
universe. These results have been verified, and made numerically more accu- 
rate, using computer simulations. However, the initial domain wall networks 
do not satisfy the static equation of motion and, therefore, must evolve in 
time. The role of dynamics in propagating or erasing these initial conditions 
is discussed in the following section. 



To investigate the evolution of the initial domain wall networks described 
above, we choose to follow Press, Ryden and Spergel (hereafter PRS)[4]. The 
dynamics of the scalar field, </>, are determined by the equation of motion 



3 Evolution of Domain Walls 



d 2 (f> 2dlnad(f> 




(24) 



drj 2 Tj d In tj dvj 
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Figure 4: The number of links across which a domain wall falls per lattice site. 
Points (a) are the 3d percolative predictions in the regime 0.175 < p < 0.25, 
with their associated uncertainty (~ 10%). Line (b) shows the measured 
values (averaged over 20 three dimensional samples) for < p < 0.5. 

where r\ is the conformal time co-ordinate, a is the scale factor of the uni- 
verse (a ~ j] in the radiation era, and a ~ rj 2 in the matter era), V is 
the scalar potential and the spatial gradients are with respect to co-moving 
co-ordinates. The scalar potential determines the topology of the vacuum 
manifold. Instead of the periodic potential given in 1, we will use the quartic 
potential 



This potential has two degenerate vacua, <f> = ±</>o, separated by a potential 
barrier Vo and closely approximates the cosine potential in the region between 
any two minima. 

Following PRS we can define a physical domain wall thickness w given 



The ratio of the wall thickness to the horizon size at the time of the ph 




(25) 



by 




(26) 
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transition 

wn 1 din a 

W - 



a{rio) i]o dlni] 



(27) 



Vo 



then sets 770, the conformal time of the phase transition and the time at which 
we begin the simulation. 

In order to avoid some technical difficulties, we will use a generalization 
of equation (24) given by 

ft + -'^ d 4 - = -«4t- (28) 

or\ l i] ami] or] o<p 

Henceforth, a = f3 = 2 which reproduces equation (24) will be replaced with 
a = 3, f3 = 0. We refer the reader to PRS for a full justification of this 
change. 

We evolve equation (28) on a regular Cartesian grid with periodic bound- 
ary conditions. Our finite difference scheme is second order accurate in both 
space and time, with the lattice equations 

1 Ari din a . 

2 i] din i] 

^ 2( t > ) . . , = <t>i+l,j,k + <^-l,fc,fc + <t>i,j+l,k + <t>i,j-l,k 

+ <f>i,j,k+i + <f>i,j,k-i - 6<f>i,j,k, (30) 

(31) 



Here, subscripts refer to the spatial lattice co-ordinate, superscripts refer to 
the time co-ordinate, and <f> = d(f>/dr]. The spatial grid size will be chosen 
to be Ax = 1. In this paper, we will set <f> = 1. The scalar field initial 
conditions are then chosen using the prescription described above for various 
bias probabilities, p. That is, in the following we will use percolation theory 
with allowed field values of ±1 at any lattice site. We will also choose the 
initial field "velocity", </>, to be zero everywhere on the lattice. 

Simulations were run in the radiation dominated epoch, with a = (77 / 770) 
and the initial time, rj = 1. We chose a wall thickness w = 5, making the 
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ratio Wo = 5. This value was used to ensure that the wall thickness was 
well above the lattice resolution scale (recall Ax = 1), while ensuring that 
for most of the dynamic range of the simulation, the wall-wall separation 
exceeded the wall thickness. The results are the following. 

We present the evolution of the energy density of the network of domain 
walls in the radiation dominated epoch. As each simulation is evolved, the 
comoving wall area is determined (according to the prescription of PRS) and 
a plot of this area, A, per co-moving volume V versus conformal time is 
produced. The simulations are run on a cubic lattice with L 3 sites where 
L = 128. They were run until a time rj = 128, or until no more walls 
remained in the box. The domain wall thickness was again set to w = 5, 
leaving us with only modest dynamical range in which to follow the network 
evolution. 

A plot of A/V for these runs is shown in Figure 5. The self-similar 
evolution seen in PRS for the p = 0.5 case is well reproduced in the time 
range 2w < rj < L/2. Fitting the scaling portion (10 < rj < 64) of the curve 
to the power law 

A/Vocf } (33) 

we find v = —0.89 ± 0.06. This is to be compared with the value found 
in PRS of v = —0.92 ± 0.06 for their three dimensional simulations in a 
matter-dominated epoch. 

Moving away from the p = 0.5 case, one sees a dramatic departure from 
self-similar scaling. In each case there is an exponential cut-off in the ratio 
A/V at some characteristic time. For the cases of p close to 1/2, that is for 
0.49 < p < 0.5, we find that the curves are well fitted by a function of the 
form 

A/V oc jfe~^. (34) 
However, for the cases where p < 0.49 a simple exponential suffices 

A/V oc e"^. (35) 

Values for v and r/, averaged over 5 runs for each value of p } are given in 
Table 1. 

For p close to 0.5 the domain wall network appears to enter a quasi-scaling 
regime in which A/V scales oc if ^ before eventually being exponentially cut- 
off at rj = fj. In particular, for the bias p = 0.499 we see that the network 
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Figure 5: Evolution of the comoving area, A, of domain walls per volume, 
V, with conformal time, 77, in three dimensional runs. For each bias, p } we 
show the evolution of one realization. 
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Table 1: Fits to the plots of A/V against rj for different initial bias probabil- 
ities, p } in three dimensions using the functional forms (34) and (35) given 
in the text. +: for the p = 0.499 case we report the best fit to equation (34), 
although the late time decay was found to be somewhat steeper than an 
exponential. 
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scales exactly as the p = 0.5 case in the epoch 2w < rj < 30, before the 
exponential decay is established. As p — > p c = 0.311, however, fj rapidly 
approaches the resolution size of the grid, and no evidence of early scaling is 
seen. This behavior continues as p drops below the critical threshold. 

To conclude, in the three-dimensional simulations we see persistent scal- 
ing behavior precisely at p = 0.5. For p below 0.5 but above p ~ 0.49, 
we see scaling for a finite time which is then exponentially cut-off at some 
conformal time fj. The value of fj } which becomes very large as p — > 0.5, 
decreases rapidly as p approaches 0.49. For p below this value no scaling 
behavior is seen and the behavior is well described by a simple exponential 
for all conformal time. 

4 Discussion 

We conclude that domain walls formed during a biased non-equilibrium phase 
transition do not, in general, persist in time. For p / 1/2 they scale for a 
relatively short time before decaying away exponentially. Only when p = 1/2 
exactly does persistent scaling behavior set in. It follows that biased domain 
walls evade the "no-go" arguments that apply to walls in thermal equilibrium. 
Since domain walls exhibit localized energy density, they will act as the seeds 
for the formation of large scale structure. A first attempt to analytically 
calculate the precise form of this structure, using percolation theory and the 
spherical matter collapse model, has been presented in [5]. The lesson learned 
there is that domain wall induced structure leaves a distinctive signature 
which could well open a cosmological window onto the microscopic structure 
of particle physics. 
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